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Abstract 

Electrostatics plays a key role in biomolecular assembly. Oppositely charged biomolecules, for 
instance, can co-assembled into functional units, such as DNA and histone proteins into nucle- 
osomes and actin-binding protein complexes into cytoskeleton components, at appropriate ionic 
conditions. These cationic-anionic co-assemblies often have surface charge heterogeneities that 
result from the delicate balance between electrostatics and packing constraints. Despite their im- 
portance, the precise role of surface charge heterogeneities in the organization of cationic-anionic 
co-assemblies is not well understood. We show here that co-assemblies with charge heterogeneities 
strongly interact through polarization of the domains. We find that this leads to symmetry break- 
ing, which is important for functional capabilities, and structural changes, which is crucial in the 
organization of co-assemblies. We determine the range and strength of the attraction as a function 
of the competition between the steric and hydrophobic constraints and electrostatic interactions. 
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I. INTRODUCTION 



The structure of biomolecular assemblies of cationic and anionic molecules containing hy- 
drophobic groups is controlled by the architecture of the molecules and the intermolecular 
interactions of the components. The co-assemblies can have spherical, cylindrical or planar 
symmetries. Oppositely charged lipids, for example, co- 

or into vesicles|2], while mixtures of cationic and anionic peptide amphiphiles capable of 
beta sheet formation co-assemble into fibers ^,0,0]. The stability of these assemblies is de- 
termined by their surface properties, which in turn is the result of the competition between 



the different type of interactions, similar to proteins j6|. The adsorption of charged lipids 
onto oppositely charged surfaces [T]] is a simple example of planar co- assembly with com- 
peting interactions where electrostatics favors homogeneous surface coverage and the lipids 
hydrophobic tails favor bilayer formation with local excess charge leading to surfaces with 
charge inhomogeneities. These charge heterogeneities are expected in many cationic-anionic 
co-assemblies including actin and microtubular associating proteins (MAP) complexes^], 
as well as in nucleosomes. Herein we analyze interactions among co-assemblies with sur- 
face charge heterogeneities. We show that co-assemblies with homogeneous surfaces interact 
weakly and induce small changes in the surface heterogeneities of neighboring assemblies 
suggesting that these co-assemblies have decreased functionality. Meanwhile, co-assemblies 
with large surface heterogeneities interact strongly and are capable of large structural mod- 
ifications. 

What controls the surface heterogeneities in cationic-anionic co-assemblies? The net in- 
compatibility among chemically different charged components, which in water can be due 
to different degrees of hydrophobicity, promotes macroscopic segregation of the different 
components. Meanwhile, electrostatics favors correlated ionic crystal structures. The com- 
petition of these interactions result in the formation of surface charge patterns ^^Is, IQ, 11| 



in aqueous interfaces due to the high permittivity of the water. In dense media charges are 
paired, while in water they can be dissociated. At the interface between a dense media and 
an aqueous solution, charged patterns of different length scales are possible because the mean 
permittivity of the two media is sufficiently high. The formation of large domains, which 
decrease the interface among the components is restricted due to the penalty associated with 
charge accumulation. Surprisingly this penalty may be overcome in the presence of other as- 
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semblies, by correlating the oppositely charged domains on the neighboring assemblies. Our 
computer simulations and theoretical analysis show that the attraction mechanism among 
co-assemblies is strongly dependent on the degree of compatibility among the co-assembled 
molecules. Surfaces of co-assembled molecules with large degrees of incompatibility have 
strong attractions which decay exponentially with separation distance, and the range and 
strength grow as the degree of incompatibility among the components increases. Compat- 
ible co-assembled components have both weak correlated interactions at short distances, 
and weak dispersion type forces at longer distances similar to the attraction among charged 



surfaces via their counterions 
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It longer 



Formation of surface patterns is not only the result of the interplay of short- and long- 
range forces. The geometry of the assembly also plays a key role. Cylindrical geometry, 
where the charged units are confined to the surface of a cylinder, is common in biomolecu- 
lar assemblies. The appearance of large domains on the surface of cylindrical co-assemblies 
induce restrictions on the size and symmetry of the domains that affect their physical proper- 
ties. Such assemblies have attracted much attention as they provide a possible way to create 
tunable functiona_l^biomolecular materials, such as stable vesicles for drug delivery and 



bio-active fibers 



Bundle formation is highly prominent in systems of fibers. Even weak dispersion forces 
lead to aggregation. In many biological systems the attractions, however, are strong and 
specific. Here we show that electrostatic interactions are important by modelling cylindri- 
cal assemblies of molecules with chemically incompatible oppositely charged head groups 
exposed onto the surface. This is a coarse grained model for cationic-anionic amphiphile 
co-assemblies or for self-attracting (hydrophobic) charged particles adsorbed on the surface 
of oppositely charged fibers such as actin-MAP complexes where the charged proteins effec- 
tively change the local properties including local charge density and degree of hydrophobic- 
ity 1^. We consider each aggregate as a stable structure, thus constraining the cyhndrical 
geometry and allowing molecules to move only on the surface. The incompatibility among 
oppositely charged co-assembled molecules per thermal energy KbT, denoted by x (expected 
to be a complex function which includes all the terms opposing mixing of molecules), and the 
electrostatic interactions determine the surface charge heterogeneities. The study of these 
interacting assemblies sheds light on the mechanism of non-specific interactions in many 
biological systems with surface compositional heterogeneities. 
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The interaction among assemblies with surface compositional heterogeneities have some 
analogies with the attraction between charged membranes or polyelectrolytes in the solution 



of neutralizing counterions 
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22]. However, there is a signifi- 



cant difference. In particular, in case of charged surfaces in the solution of neutralizing 
counterions, the size of the Wigner crystal cell is determined by the conditions of surface 
electro-neutrality and A oc \/QI {ea\ where ea is a surface charge density and Q is the 
charge of condensed counterions 16|. The fluctuations of condensed ions and structural 
correlations induce electrostatic attraction. On the other hand, in case of the surfaces with 
compositional heterogeneities, the size of the Wigner crystal cell is determined by a balance 
of electrostatic and interfacial energies and A oc a/x7^- should be also noticed that 
the area of the domain determines the effective charge Q oc L^, while charge of condensed 
counterion remains the same. 

In this paper we study the interaction among cylindrical assemblies with surface com- 
positional heterogenous. We show that the compositional heterogeneities of two interacting 
cylindrical assemblies are correlated and the size of heterogeneities grows as the assemblies 
approach each other. This leads to a strong attraction. We show that the effective potential 
of interaction depends on the degree of surface compositional segregation and the separation 
distance among assemblies. The range and the strength of the attraction increase as the 
degree of the net incompatibility of the cationic-anionic components increases. 



II. SIMULATION METHOD 



We consider each aggregate as a stable structure, thus constraining the cylindrical ge- 
ometry and allowing molecules to move only on the surface. Within the framework of the 
primitive model, each aggregates is composed of a stoichiometric mixture of equal number 
= of positively and negatively charged units of equal absolute charge \Q+\ = \Q~ \ = 1 
and diameter a = 1. The cylindrical assemblies are placed into the box x L^^ x in the 
middle of the FZ— plane along 2;— axis, where Ly = Lx + 2Rc + D. We study the interaction 
of two parallel cylindrical assemblies of radius Rc and length separated by the distance 
2Rc + D, where D gives the distance between ions on the faced surfaces of the neighbor- 
ing assemblies. From here, due to the fixed radius Rc, we will use only D to indicate the 
separation distance. 
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The surface of each cyhnder is filled with units in such a way that all units are placed into 
the knots of a triangular lattice of period a (Fig. and thus the number surface density 
IS cr = 2/V3a2. The short -ranged interaction between units is taken into account via a 
simplistic discrete van der Waals potential, due to the discreet model. The net degree of 
compatibility among different components, i.e. positively and negatively charged units, is 
controlled by Flory-Huggins parameter x = l^ij ~ ^{^a + ^jj)]/kBT, where Eij represents 
a pair interaction energy among i and j ions in fc^T units. Compositional incompatibility 
results in the formation of segregated domains, whose growth is inhibited by electrostatics. 
We use Debye-Hueckel potential to describe the electrostatic interactions 

u{r)/kBT = e—- (1) 

where i = /ATreoksT is Bjerrum length, 60 is the average dielectric constant of the media 
in units e'^/Aira and n = AniY^Cizf is the inverse Debye screening length, q is the concen- 
tration and Zi is the valence of the i—th neutralizing counterion in the solution. Taking into 
account a biological origin of the system, i.e. the large head groups size (~ lOA) and high 
average dielectric permittivity of the aqueous solution (~ 80), the Bjerrum length i = 0.2a 
is considered. Due to the assemblies are electroneutral, counterion condensation in not 
expected in the wide range of temperatures and, thus, we choose Debye length a/hz = 25. 

We report standard canonical Monte Carlo simulations following the Metropolis scheme 
for various values of x G [0.5, 12], Rc/a = 2 — 5, = and L^/a = 100. Simple moves in 
the phase space are performed by exchange of two randomly chosen particles. Each system 
is equilibrated during 10^ MC steps per particle and another 10^ MC steps are used to 
perform measurements. The equilibration process is accompanied by a gradual decrease of 
temperature (temperature annealing) from T^ax = 10 to T^m = 1- 

III. RESULTS 

The competition between the net short range repulsion among dissimilar components and 
the electrostatic penalty associated with the charge accumulation results in the formation 
of the surface composition heterogeneities. With increase in x domains size grows, however 
the electrostatics restrict the possibility of macroscopic segregation. The steric commen- 
surability between the surface charge patterns applies geometrical restrictions on the size 
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and the symmetry of the domains, and together with the long-range electrostatic inter- 
actions induces periodic arrangement of domains of opposite charge, which stabilizes the 
co-assembled structure in analogy with the Wigner crystals of ions 23[. Figure 121(a) shows 
that with increase in x the normalized energy of the system VjxN is decreasing for both, 
nearby [D jo = \) and separated assemblies. The energy of the nearby assemblies is notably 
smaller, suggesting a strong attraction among aggregates. 

With increase in x the system undergos transition from the isotropic to segregated phase 
that is indicated as a broad peak on the heat capacity (Fig. |21 (b)) at x denoted by X\- 
For nearby assemblies the peak is higher and it is shifted to smaller x values, suggesting a 
promotion of the compositional segregation by the assemblies cross-interaction. 

A. Theoretical model 

Due to complexity of analyzes of cylindrical assemblies, we start by considering two 
interacting planar surfaces. In the limit of small net incompatibility, the contribution of 
density fluctuation to the mean free energy can be analyze by means of Random Phase 



Approximation (RPA) 2^ . In the limit of small density fluctuations around an homogenous 
state of mean density p equal to the number fraction of positively and negatively charged 
components p = p+ = p_. The total free energy F including one loop corrections (OLC) 
can be obtained from the partition function Z = e"^/'^^-^. The partition function of the 
system in terms of the Fourier components of the density Pg reads 




where = e~^°/^^^ is zero mode partition function (g = 0) with no fluctuation contribution. 
Pi ~ (Pq^Pq)^ where and p^ are the Fourier components of the composition fluctuations 
around the mean value, Ap*(r) = p*(r) — p, of cylinder i = 1,2, respectively. A is the inverse 
of the density correlation matrix S with elements Sij{q) = (^p^^p^^^. The diagonal elements 
aii{q) of the matrix A can be estimated according to Cahn-Hilliard free energy of a neutral 
system j^l with additional gradient (interfacial energy) and electrostatic terms 

1 { Ci] ^ 

= p(l-p) ~^^ + ^^ +^<yuu{q) (3) 



6 



where the electrostatic potential 

uu{q) = i / —e-'^dr = -^===. (4) 

On the other hand, non-diagonal terms, aij{q, D) = AaUij{q, D), include only the energy of 
the electrostatic interaction between the separated surfaces 

> 

27r£ 



UiJq,D) =i / e'^^dr 



e 



where D is the distance between the surfaces. 



B. Domains polarization 

In the limit D ^ oo, when assemblies do not interact, the competition between elec- 
trostatics, AaUii{q), and interfacial, x(9'^)^/2, energies result in the formation of favor- 
able density fluctuations of finite wave length g* and appearance of a peak in Sii(g) at 



q* ~ (87r£cr/xa^)^/^(Fig. 121 (a)). Besides, detailed analysis o: 



tures restricted on the surface of cylinders 



and planes [h] reveal many interesting final 



binary cationic-anionic mix- 



temperature effects as well as stripe structures along the assembly at lower temperatures. 

As one cylinder approaches another, long range electrostatics correlates surface compo- 
sitional heterogeneities on the neighboring assemblies to minimize the electrostatic energy. 
In closed vicinity, D < R^, domains of opposite charge are located in front of each other 
and electrostatic energy of the cross interaction effectively compensates penalty associated 
with charge accumulation. That results in the polarization or enlargement of the domains. 
The peak position is shifting in direction of small q values as D decreases (Fig. El (a))- In 
the limit q*D > 1, the dependence of the domain size on the distance D among assemblies, 
calculated form the dsu/dq = 0, has simple approximation 

g*^g:(l-g:De-^S^)^/l (6) 

At the same time, the magnitude of the peak S{q*) is growing with decrease in distance 
that indicates strengthening of the compositional segregation. Figure El (b) shows the OLC 
critical conditions, obtained from the numerical solution of equation = 0. 
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To analyze the domains polarization and compare results of the computer simulation and 
theory we calculate the static structure factor along the assembhes faced surface 




where A^;^ = Zi — Zj for all ion pairs with equal X and Y components, Xi = Xj and 
Vi = Vj, where = {xi,yi,Zi) is a cartesian vector. The competition between electrostatic 
repulsion and interfacial energies results in the formation of most favorable fluctuations of 
finite wave length A = 27r/g* and the appearance of a peak in Sz{q) at q*. Figure |3](b) shows 
the dependence of q* versus D for different values of x- With decrease in the separation 
distance the q* is also decreasing, indicating the domains enlargement. For x < X^^\ the 
q* always has a final value, while for x > X^^^ fo^' nearby assemblies the domains size 
becomes in order of the assemblies length (Fig. 0] (a)) indicating a very strong segregation. 
It is important to underline, that in this case cylindrical geometry results in breaking of 
symmetry of surface domains, due to A ~ R^, and then domains grow only along the 
assembly. Polarization of domains is not only a new and interesting phenomenon that may 
be used to explain aggregation of different assemblies or biomolecules into bundles. It also 
affects the effective potential of attraction, due to the dependence of domain size A = 2n/q* 
on the separation distance D. Polarization of domains shows another significant difference 
in the nature of interaction among assemblies with surface compositional heterogeneities 
and charged membranes or polyelectrolytes in the solution of neutralizing counterions. 



C. Electrostatic attraction 



The correlation and polarization of surface domains decreases the energy of the system, 
suggesting the attraction among the assemblies. The strength of the attraction can be 
calculated from the dependence of energy U / xN on the distance D among assemblies, as 
AU{D)/N = {U{D) — Uo)/N, where Uo is the energy of two separated assemblies. At the 
same time, the energy of two interacting planar assemblies can be analyzed analytically as 



^,d{F/T) 



^E,, = -T'^^, (8) 

where F is the OLC free energy (Eq. E)). Figures |S1 (a) and (b) show both of the energies for 
different values of the net incompatibility x- 

8 



In the limit of small net incompatibility, x < Xo, when the surfaces are nearly homo- 
geneous and isotropic, the attraction is originates by cooperative correlations of charge 
fluctuations leading to a long-ranged power law attraction, [22]. With decrease in the 
distance the density fluctuations are suppressed due to the lateral repulsions among 
ordered domains. That leads to stronger attraction. On the other hand, with increase in 
the degree of net incompatibility x the strength of the attraction is also growing, due to 
the increasing magnitude of the density fluctuation. However, as we will show later, in the 
limit of strong incompatibihty, x > Xo, we find an unexpected slowly decaying exponential 
attraction. 

In the limit of strong net incompatibility, when surface domains are well pronounced, 
charged domains correlate themselves to minimize the energy and form one dimensional pe- 
riodic array. Therefore, the ground state configuration is a pair of two cylindrical assemblies 
with a charged surface pattern of positive and negative domains in one assembly in front of 
an oppositely charged domains in the second assembly. In that case, the attraction energy 
can be approximated as 



This approximation should work in the limit of thin cylinders, when domain size A is large 
then the radius of the assembly Rc. On the other hand, when A < Rc, cylindrical assemblies 
can not be considered as one dimensional objects and a planar surface represents another 
limit. In this case, the attraction energy scales as 



It should be noted here that the domain size A is expected to depend on the distance between 
assemblies. 

A possible fitting of the results (Fig. El (a)), which sheds some light on the origin of the 
attraction and also includes the dependence is 








(10) 



AU^':>/N = U^,^/N-P,D 



(11) 
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where d is the space dimensionahty, U2d/N = —Eie~^^^ and Uid/N = —Eie~^^^ /\/D, 
and Ei,E2, Pi and P2 are fitting parameters. Taking into account the D dependence of the 
domain size A, we can expect more complex dependence on the separation distance. 

Figures El (c-f) shows the dependencies of all fitting parameters on x and compares the 
corresponding fitting parameters with the parameters we computed analytically. When 
the system is nearly isotropic, the attraction always includes a short-range (exponential) 
component and a long-range (power law) component. The short-range term {Ei, E2) is found 
in the whole range of x, suggesting correlations even for small degrees of net incompatibility. 
Meanwhile the long-range term (Pi, P2) is found only for x < X^^\ when the system is mostly 
isotropic. At the point x^^^ the system undergoes a crossover transition from isotropic to 
periodic phase and that affects the nature of interaction of the aggregates. For x < X^^^ both 
parameters Ei and Pi, those determine the strength of each component of the attraction, 
exponentially grow with x, while for x > X^^\ the long-range term diminishes and the short 
range term changes the dependence on x- In the strong segregation limit, x > Xo, the long- 
wavelength density fluctuations and therefore the power law component of the attraction 
are suppressed by strong domain correlations and the x~dependence of the exponential 
term changes with respect to the x < Xo limit. The multiplicative coefficient Ei remains 
nearly constant, while the exponential decay E2 starts decreasing suggesting a growth in 
the range of attraction. This is due to the fact that the domain size on the face of the 
interacting cylinders also grows (Figs. E] (a)). In the limit of large x > X^^\ since there no 
long-ranged order in ID, the long-wavelength density fluctuations are suppressed together 
with the long-range component of the attraction. 



IV. CONCLUSIONS 



Compositional heterogeneities and the changes induced by nearby assemblies via cor- 
relations and polarization, are key in the design of functional supramolecular assemblies. 
Correlated domains effectively cancel the electrostatic penalty of each other and amplify 
compositional segregation. This leads to strong attractions. The cylindrical geometry re- 
sults in the breaking of the symmetry of the surface patterns, which is important in triggering 
functional properties. For example, in systems with competing interactions it can explain 
aggregation of co-assemblies into clusters of different symmetries. The bare length scale 
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of the charge surface heterogeneity and the degree of polarization determine the symmetry 
of the bundles of fibers perpendicular to their long axis and unfolds the origin of hexag- 
onal, square or linear (ribbons of fibers) symmetries of supramolecular structures. These 
effects are of general importance in the organization of biological assemblies. Our results not 
only explain the mechanism of the attraction between both planar and curved surfaces with 
charged heterogeneities, but also give insight to the understanding of electrostatic effects in 
molecular biology 0, 1^6 1 . 

A recent experiment [7| shows the strongly correlated attraction mechanism discussed 
here among negatively charged surfaces (mica) coated with cationic hydrophobic surfactants 
(DODA). The surfactants segregate into charged bilayer domains. This leads to a strong 
attraction among coated parallel surfaces due to the growth of the bilayer charged domains 
and the correlations between them and the bare mica domains on the neighboring surface. 

In summary, we have studied interaction among assemblies with surface compositional 
heterogeneities. The transition from isotropic to ordered phase leads to formation of charged 
domains along the cylindrical assembly. The attraction appears as a result of correlations and 
polarization of charged domains. We show that the interaction among assemblies belongs 
to the general class of fluctuation-induced forces and includes long- and short-range term. 
Formation of periodic domains along the cylinder leads to the long-range order of charged 
domains that suppresses long-wavelength density fluctuations and, thus, long-range com- 
ponent of the attraction. Attraction correlates and polarizes charged domains. Our results 
reveal the importance of charge heterogeneities in cationic-anionic co-assemblies of complex 
molecules and suggest a promising strategy for fabrication of assemblies with predictable 
surface charged patterns for developing functional biomolecular assemblies. 

This work is supported by NSF grant numbers DMR-0414446 and DMR-0076097. 
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FIG. 1: Snapshots of typical configurations of (a) isotropic phase {x = 1-5) and (b) striped state 
with defects (x = 10.5). D/a = 3, Rc/a = 5 and Lja = 100. 



FIG. 2: The dependence of (a) normahzed energy U /xN and (b) dimensionless, normahzed by N, 
heat capacity versus x fo^' different separation distances D. 

FIG. 4: (a) Snapshots of typical configurations of surface domains on the faced sites of separated 
and nearby assemblies for x = 3.0, 6.0 and 9.0. (b, c) Dependence of q* on the separation distances 
D calculated in the simulation and numerically {dsn/dq = 0). 



FIG. 5: (a) Normalized energy of interaction (points) for two cylindrical assemblies separated by 
the distance D for different degrees of x- Solid lines present results of fitting of normalized energy of 
interaction (Ea. [TT|) . (b) energy of interaction calculated according to Eq.|Hl (c-f) fitting parameters 
Ei,E2,Pi and P2 versus x- 



FIG. 3: (a) The density correlation function sii{q) for different separation distances D and (b) the 
RPA critical function of D. 
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